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ABSTRACT 

In  this  paper,  we  introduce  a new  continuous  method  for  the  general  convex  problem  (1).  Our  method  overcomes 
all  the  shortfalls  of  the  previous  methods,  yet  maintains  strong  convergence.  As  pointed  out  in  [2]  that  the  projection  model 
(3)  is  difficult  or  impossible  to  be  implemented  for  a general  convex  constraint,  our  approach  first  transforms  problem  (1) 
into  a monotone  variational  inequality  problem.  Then,  a continuous  method  is  introduced  for  this  variational  inequality 
problem.  Our  continuous  model  differs  from  the  one  in  [3]  in  two  aspects.  First,  we  employ  a merit  function  in  our 
continuous  model.  The  structure  of  our  model  is  similar  to  the  one  in  [23,  9].  Second,  the  bound  for  the  scaling  parameter 
(called  do  in  [3])  is  constant  in  [3],  but  it  is  variable  in  our  continuous  model;  in  particular,  an  ODE  for  this  variable  as 
function  of  time  is  established. 

KEYWORDS:  Convex  Programming,  Variational  Inequality,  Continuous  Method 

1.  INTRODUCTION 

In  this  paper,  we  consider  the  following  convex  programming  problem: 

(CP)  min  eS„  / (x)  (la) 

Such  that  g(x)  < 0 , g : R"  — > R'"  (lb) 

Where  f(x)  and  g(x)  are  convex  functions  and  have  continuous  first-order  derivatives.  Obviously,  for  convex 
problem  (1),  every  local  solution  is  also  a global  solution.  For  problem  (1),  we  define 

n0={*€r|g(*)<o} 

as  the  feasible  set  and 

={xe  R"  |jc  is  an  optimal  solution  of  (1)} 

as  the  solution  set.  Throughout  this  chapter,  we  assume  that  there  exists  a 

finite  A*  £ T2()  and  that  the  Slater  condition  [i.e.,  3 an  x such  that  g(x)  < 0,  [5]]  is  satisfied.  Thus,  the  KKT 
conditions  for  problem  (1)  can  be  written  as 

Vf(x)  + (g'(x))T=  0,  (2a) 

T - 0 , g(x)  < 0 , yrg(x)  = 0 (2b) 
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Even  though  continuous  methods  can  be  traced  back  to  the  1950s  [1],  the  most  significant  theoretical  contribution 
to  continuous  methods  comes  from  Hopfield  [16-18].  In  recent  years,  the  restriction  on  the  merit  function  being  a 
Lyapunov  function  has  been  removed  to  allow  wider  applications  [25,  20,  10]. 

As  in  Flam  [7]  the  ODE 

dx(t)/dt  =-lx(t),  X >0. 

The  solution  of  the  above  ODE  is 

x(t)  = x0e 

Thus,  it  is  easy  to  see  that  x(t)  — > Oas  t — >co.  But  x(t)  will  never  reach  x*  = 0 in  any  finite  time  if  x,  ^ 0 . 

Therefore,  the  conditions  and  result  of  theorem  3.2  in  [7]  do  not  hold.  The  first  convergence  proof  of  a continuous  model 
for  the  general  convex  problem  (1)  was  given  in  [2].  The  ODE  suggested  in  [2]  is  based  on  the  gradient  projection  method 
and  has  the  following  form: 

dx(t)/ dt  = -[x-Pn  (*-V/ (.*))]  (3) 


Where  the  feasible  is  set  of  (1)  and  P 0 (•)  is  the  projection  onto  £2() . Even  though  this  method  is 
theoretically  sound,  it  has  serious  practical  difficulties,  since  the  projection  P ()  (•)  is  difficult  to  compute  a general  convex 
set  Q.0  . 


2.  PROBLEM  EQUIVALENT  TO  VI 

First,  let  us  look  at  the  following  lemma. 

Lemma  2.1.  .V'  €E  f2(|  <=>  3 ay  G R'"  such  that  x and  y satisfy  (2). 

Proof.  By  [22]  for  =>  and  Theorem  in  [8]  for  <=. 

Now,  let  us  define 

' N rVf(x)  + (g'(x))Ty^ 

s k~S(x)  )’ 

and  consider  the  following  VI  problem: 

(VI (f 1 , F))  Find  U*  e Q such  that  (u  - U*)T  F(u *)  > 0 , Vw  € Q 


Where 


Q = < U 


yyJ 


XG  R",vg  R”  [ ,/?:"  ={_vg  R:  y > ()} . 


Then,  the  following  lemma  shows  that  F (u)  defined  in  (4)  is  monotone. 


(4) 


(5) 
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Lemma  2.2.  F (u)  in  (4)  is  monotone  on  Cl . 

Proof.  For  any  U = ( XT  ,yT)  , U = (x T , y T )T  € Cl  we  have 


(u  - u)T (F (u)  - F (u)) 


fx-x^T 


yy-y) 


(Yf(x)  - V/(x)  + (^^))r  y - ( g'(x))T  y) 


= (x-x)T (Vf(x)-Vf(x))  + yT {g(x)~ g(x)~ g'(x)(x  -x)}  +yT {g(x)~  g(x)~  g'(x)(x-x)} 

(6)  Since  f(x)  is  convex,  we  have 

(x-x)T(Vf(x)-Vf(x))>  0 (x-x)>0  (7) 


Since  g(x)  is  convex,  from  Proposition  4 in  [21],  we  have 

{<?(*) + ■?'(*)(*--*)},  (8a) 

g(x)>  {^(a)  + g\x)(x-x)}  (8b) 


Replacing  (7)  and  (8)  into  (6)  and  noticing  that  y > Oand  V > 0 , we  have 


(u-u)t(F{u)-F{u))>  0. 

This  completes  the  proof  of  Lemma  2.2. 

Lemma  2.2  indicates  that  (5)  is  a monotone  VI  problem.  Let  Cl*  be  the  optimal  solution  set  of  (5). 
Then,  we  have  the  following  lemma. 

Lemma  2.3.  u satisfies  (2)  if  and  only  if  U E Cl  * . 

Proof.  Let  U = ( XT  ,yT)T 


(=>)  Since  u satisfies  (2),  then  U < 


Cl  and,  Vm  = (xT  ,yT)T  1 


Cl,  we  have  (u  - u)T  F (u)  = -(y  - y)T  g (x) 


= -yTg(* ) (9) 

Obviously,  y > 0 .Then,  g(x)  < Oand  (9)  imply  {u  —u)TF(u)  > 0,  Vw  € Cl . 

(<=)  Since  we  Cl  * , wehave  y > 0 and 


(A-A)r[V/(A)  + (g/(^))r  >’J  -(T-T)r<?(^)^0  Vw  e Cl 


(10) 


Let  the  vectors  y — y and  X in  R"  . Then,  (10)  implies  that 

vf{x)+(g'{x)fy  = o 

Consequently, 
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(ID 


y = V + 6- , 1=1,2, Ill  where  ej  is  the  ith  column  vector  of  the  identity  matrix  Im.  Then,  (11) 

becomes  —ejg(x)  > 0 , i =1, 2, m . 

T 

These  inequalities  imply  that  g(x)  < 0.  Now,  we  need  only  to  verify  that  y g(x)  = 0.  Since  y > 0 and  (11)  is  true 
for  all  y > 0 , setting  y = 0, 2 y respectively,  we  obtain 

T T 

y g(x)  > 0 and  - y g(x)  > 0. 

T 

This  implies  that  y g(x)  = 0.  Therefore,  u satisfies  (2).  This  proves  Lemma  2.3. 

3.  PROJECTION 

For  any  V € Rm ' " , we  denote  by  Pa  (v)  the  projection  of  v on  £2  under  the  Euclidean  norm.  In  other  words, 
pa(v)  = argmin|||v  — m||m£  £2}  is  the  unique  solution  of  the  minimization  problem  minjjv  — w|| \u  € £2}.  Notice 
that  the  projection  on  D.  [5]  formulated  from  the  convex  programming  problem  is  trivial. 

A basic  property  of  the  projection  mapping  on  a closed  convex  set  is  that 

(v  - Pn  (v))T  (u  - Pn  0))  < 0 , Vv  e /?" , Vm  e Q (12) 

Since  the  early  work  of  Eaves  [6]  and  others,  we  know  that  a variational  inequality  VI  {Q.,F)  is  equivalent  to 
the  following  projection  equation: 

u=Pn[u-F(u )]  (13) 

In  other  words,  solving  VI  (Q.,F)  is  equivalent  to  solving  the  nonsmooth  equation  (13),  i.e.,  finding  the  zero 
point  of  the  residue  function 

e(u)  = u - Pa[u  -F(u)]  (14) 

Because  e(u)  is  continuous  and  since 

e(u)  = 0 <=>  Vwe  £2* 

Let  M*g£2*  be  a solution.  For  any  U £ R" , Pn [ it  — F(u)\&  Q. . It  follows  from  (5)  that 

F(u*)r{Pa[u-F(u)]-u*}>0  , Vwe  R”  (15) 

Setting  V — U — F(u)  and  U — U * in  (12),  and  using  the  notation  e(u) , we  obtain 

{e(u)-F(u)}T{Pn[u-F(u)]-u*}>0  , Vwe  Rn  (16) 
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Provided  that  the  mapping  F is  monotone  with  respect  to  Q.  , we  have 

{F(Pa[u-F(u)])-F(u*)}T{Pa[u-F(u)\-u*}>  0,  Vue  R"  (17) 

(15)  Follows  from  the  definition  of  variational  inequalities;  (16)  follows  from  a basic  property  of  the  projection 
mapping;  (17)  follows  from  the  assumption  of  monotonicity  of  the  mapping  F. 

Besides  e(ll)  [7], 

d(u)  = e(u)-{F(u)~  F(Pn[u  - F(u)])}  (18) 

Adding  (15),  (16),  (17)  and  using  the  notation  e(u)  and  d(u  ) we  have 

d (u)T  ((u  - u*)  - e(u))  > 0 , Vu  e R" 
and  it  follows  that 

(u  -u*)T  d(u)  > e(u)Td(u),Vu  e R"  (19) 

Lemma  3.1.  For  Le  (0,1)  , if  the  following  condition  holds: 

\\F(u)-F(Pn[u-F{u)])\\<L\\e(u)\\,Vue  Rn+m  (20) 

then  we  have 

e(u)T d(u)  > (1-L)||e(w)||2  (21) 

Proof.  From  (18)  and  (20),  it  is  easy  to  see  that 

e(u)T d{u)  = \e{u)f  -e(u)T {F(u)~ F(Pn[u  - F(u)])} 

> \\e(u)\f  - \\e(u)\\ ||F(m)  - F(Pn[u  - F(«)])|| 

>(l-L)\\e(u)f 

This  proves  (21)  and  Lemma  3.1. 

For  some  positive  scalar  /?  we  consider  the  mapping  fd F and  denote 

e(u,/3)  = u - Pn[u  - /3 F (u)] 

d(u,  J3)  = e(u,  jB)  - (J3F(u)  - j3F(Pn[u  - / 3F(u)] )) . 

It  is  trivial  to  see  that,  for  any  J3  > 0 , 
e(u,/3)  = 0 <=>  e(u)  = 0 o u e Q.  * . 
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Now,  we  explore  some  properties  for  e(u,  j3)  . 

Lemma  3.2.  For  all  U G R"  and  j3  > j3  > 0 , it  holds  that 
\e(u,]3)\>\e{u,  fi)\ 

\\e(u,/3)\\/  /3<\\e(u,/3)\\/  J3 


Proof.  If  e(u,  /3)  = 0 , then  e(u.  j3)  = 0 . Therefore  the  results  hold.  Now,  we  assume  e(u,  j3)  ^ 0 . 


Let!  = 


e(u,j8) 


We  need  to  prove  only  that  1 < t < — . Note  that  its  equivalent  expression  is 
From  (12),  we  have 

(v  - Pn  (v))r  (Pn  (v)  - w)  > 0 , Vtv  G Q. 

Substituting 

w = Pn[u  - j3F(u)\ , v = u — j3F(u) 

in  (25)  and  using 

Pn[u  - J3F(u)\  - Pn[u  - ]3F(u)\  = e(u,]3)  - e(u,j3) , 

We  get 

{e(u,j3)-j3F(u)}T{e(u,P)-e(u,j3)}>  0 

Similarly,  we  have 

{j3F(u)-e(u,j3)}T  {e(u,j3)-e(u,j3)}>  0 

Multiplying  (26)  by  /3  , (27)  by  /? , and  adding  the  resulting  expressions,  we  get 
{fie(u,fi)-fie(u,/3)}T{e(u,/3)-e(u,j3)}>  0 

Consequently, 

fi\\e(u,fi)f  +y?|e(M,y5)||2  < ( j3  + jB)e(u,j3)Te(u,jB ) 


(22) 

(23) 


(24) 


(25) 


(26) 


(27) 


(28) 


(29) 
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Dividing  (29)  by||e(w,/?)||  , we  obtain 

P + Pt2  <(p+p)t 

thus,  (24)  holds  and  Lemma  3.2  is  proved. 

Lemma  3.3.  For  any  U £ ^"+"'  with  e(u)  i1  0 and  L £ (0, 1)  , there  always  exists  a small  P > 0 such  that 
\\PF(u)]-PF(Pn[u-pF(u)])\\<L\\e(u,P)l  V/?£  (Oj)  (30) 

Proof.  We  divide  the  proof  into  two  cases. 

Case  1 : U i.  Q. . Because  is  a closed  convex  set,  there  exist  an  £ > 0 and  a 8 > 0 such  that 

\\e(u,  fi)\\>  £ , V^£  (0,£)  (31) 

Then,  from  the  continuity  of  F , (30)  can  be  ensured. 

Case  2:  U £ Q Since  e(u)  ^ 0 , there  must  exist  a [3  > 0 such  that  ||e(w, /?)||  > 0 . Therefore, 
j8  >0 . 

However  U £ D.  implies 

lim^+0  \\F(u)]  - F(Pn[u-j3F(u)])  ||  = 0 

Together  with 

||e(M,/0||//?  >0 

and  (23),  this  indicates  that  (30)  is  true.  Combining  the  two  cases,  there  exists  a > 0 such  that  (30)  holds.  This 
proves  Lemma  3.3. 

The  result  of  Lemma  3.3  ensures  that,  by  reducing  /?  small  enough,  (30)  can  be  always  satisfied.  Then,  similarly 
to  (19)  and  (21),  we  have 

(u  -u*)T d(u,  ft)  > e(u, P)T d(u, P)  > (1-L)||e(w,/?)||-  (32) 

4.  MONOTONE  VI  PROBLEM 

In  the  continuous  method,  we  adopt  a technique  to  adjust  the  parameter  p so  that  (30)  and  furthermore  (32)  can 
be  satisfied.  We  define 

P = ea  (33) 

and  let 
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W ■ 


dw 

dt 


' du.  N 

dt 
da 
v dt  ) 


Let 


7J(a,  F,«r)  = w -Pa[w  -e“F(M)]  , 

r2 (u,F,  a)  = || eaF(u)  - eaF(Pn[u  - eaF(u)])f 

Then,  we  define 

du  _\-d(u,ea),  if  T2  < 0.97| 
dt  |().  otherwise 

da 

— = -rj(u,F,a ) 
dt 

Where 


ti(u,F,a ) 


Jo,ifr2  <0.97; 

[l,  otherwise 


Similarly,  we  can  define 


£{u,F,a) 


Ji,if  r2  <0.97; 

[0,  otherwise 


Therefore,  we  have 

t/vv _ (%(u,F,a)d(u,ea)' 
dt  [rj(u,F,a)  y 


In  (34),  the  right-hand-side  term  is  not  continuous.  To  make  this  term  continuous,  we  can  define 


dw  _ (£(u,F,a)d(u,ea)s 
dt  { ij(u,  F,a)  y 

Where 


Z(u,F,a) 


1,  if  T0  < 0.8 

< (0.9-r0)10,  if  0.8  <T0<  0.9 
0,  if  r0  > 0.9 


(34) 


(35) 


(36a) 
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r/{u,F,a) 


0,  ifro<0.8 

< (T0  -0.8)10,  if  0.8  < T0  < 0.9 

1,  if  T0  > 0.9 


T0(u,F,a ) 


T2(u,F,a) 
Tx(u,F,a ) 


Obviously,  the  right-hand-side  term  in  (35)  is  continuous  now. 

The  merit  function  is  given  as 

E[w)  = (1/2)||m  — M*||  , Vlt  G i?"+m  ^ere  ll* E and  is  finite. 

5.  CONVERGENCE  PROPERTIES 


(36b) 


(36c) 


(37) 


Here  we  will  discuss  the  convergence  properties  for  (35),  we  use  T()  in  place  of  T0(u,  F , CX) . 

Theorem  5.1.  For  any  W0  E R , there  exists  a solution  w(t ) of  (35),  with  tv(/  = 0)  = tv'0  and  tv(  / ) defined 
in[0,°o)  . 

Proof.  From  our  assumption  that  there  exists  a finite  X * € fi(| , we 
Know  that  there  exists  a finite  u*e  a*.  Then, 

~(u -u*)T d(u,ea),  if  T0  <0.8 

-(9  -10ro)(w  -u*)T d(u,ea), if  T0  < 0.8  (38) 

0,  if  T0  > 0.9 

For  all  cases  of  (38),  from  (32)  with  L =0.9,  it  is  easy  to  see 

d llzr  — u *||2 

■ ^<0  (39) 

dt 

(39)  Indicates  that 

u(t)e  B(u0,u*)  = {«e  R"+m  |||w  — u *||  < ||w0  —u  *||} 


d ||w  — u *|| 
dt 


dcx 

The  set  B(u0,U*)  is  a closed  bounded  set.  In  addition,  it  is  easy  to  see  that <0.  From  the  continuity 

dt 


of  d(u,ea)  , it  is  easy  to  verify  that  the  right-hand-side  function  of  (35)  is  bounded.  From  the  Cauchy-Peano  theorem,  the 
result  is  obtained.  This  proves  Theorem  5.1. 

Lemma  5.1  For  any  W0  G jet  vv(f)  be  a solution  of  (35)  with  w(t  = 0)  = W0  . Then  there  exists  a 
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t > 0 such  that 

7;  <0.8  if  t>T  (40) 

Proof.  For  some  l > 0 , if  71  = 0 , this  ensures  thatw(f)  G £2  * . Therefore,  (40)  is  true.  Now,  we  assume 
that u(t)£  n*,  Vf  > 0 . 

da 

First,  from  (35),  it  is  easy  to  see  that < 0 . Therefore,  (X  is  always  monotonically  non-increasing  in  t.  From 

dt 

Lemma  3.3  and  the  definitions  of  TJ  (u,  F , CX)  and  /3  ,we  have  that,  \/u(t)  £ £1  * , 3 CX(il)  such  that 
T0<O.S,\/a<a(u)  (41) 


From  the  proof  of  Theorem  5.1,  we  know  that  w(l)G  B(u0,U*)  , which  is  a closed  bounded  set.  This  proves 
Lemma  5.1. 

Theorem  5.2.  For  any  W()  G /?"+”'+1 , let  w(t ) he  a solution  of  (35)  with  w(t  = 0)  = R;() . Then, 
lim,^  u(t)  exists  and  lim,^  u(t)  e £1  * . 


Proof.  From  Lemma  5.1,  there  exists  a t > 0 such  that  T{]  < 0.8 , Vf  > t . From  (35),  (38),  and  (32),  we  can  see 
easily  that 


d u-u* 

dt 


< 


(l-L)\\e(u,ea)f  yt>t 


(42) 


Where  U*  G £1*  is  finite  and  L is  a constant  in  (0,  1). 

From  the  discussion  in  the  proof  of  Lemma  5.1,  we  know  that  a is  monotonically  non-increasing  and  that 


^=o.v,> 

dt 


t . 


Therefore,  1 i m ( ^c:  Coexists.  Let  a* be  the  limit.  Then,  we  define 


w*  = 


Ka*j 


From  (42),  we  have 


d w-w* 
dt 


< 


-(l-L)\\e(u,ea)f  yt>t 


From  [24]  and  (43)  we  know  that 

lim^oo  e(u(0,  ea*)=  0 or  lim,^  e(u(t))=  0 


(43) 
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From  the  proof  of  Theorem  5.1,  we  know  that  u(t)  £ B(u0,U *)  , which  is  a compact  set.  Therefore,  there  exists  a 
sequence  tk  , with  tk  — > °°  as  k — > °°  , such  that  limA  oo  U (t,  ) exists.  Let  ll  be  the  limit.  From  (42),  it  is  easy  to  see  that 

e(u,ea*)  = 0.  But e(u,jB)  = 0«  e(u)  = 0 <=>  if  £ £2*  . Therefore,  u e £1* 

By  replacing  U*  by  U in  (42),  we  have 


7.  CONCLUSIONS 

In  this  paper,  we  have  discussed  a continuous  method  for  convex  programming  (CP)  problems,  which  includes 


both  a merit  function  and  an  ordinary  differential  equation  (ODE),  for  the  resulting  variational  inequality  problem.  Our 
approach  is  different  from  existing  continuous  methods  in  the  literature.  Strong  convergence  results  of  our  continuous 
method  are  obtained  under  very  mild  assumptions.  No  Lipschitz  condition  is  required  in  our  convergence  results. 
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This  completes  the  proof  of  Theorem  5.2. 


www.tjprc.org 


editor@tjprc.org 


154 


Prashant  Chauhan 


10.  Han,  Q.  M.,  Liao,  L.  Z.,  Qi,  H.  D.,  and  Qi,L.Q.,  Stability  Analysis  of  Gradient-Based  Neural  Networks  for 
Optimization  Problems,  Journal  of  Global  Optimization,  Vol.  19,  pp.  363-381,  2001. 

11.  He,  B.S.,  Liao,  L.  Z.,  Han,  D.  R.,  and  Yang,  H.,  A New  Inexact  Alternating  Directions  Method  for  Monotone 
Variational  Inequalities,  Journal  of  Optimization  Theory  and  Applications,  Vol.  92,  pp.  103-118,  2002. 

12.  He,  B.  S.,  and  Liao,  L.  Z.,  Improvements  of  Some  Projection  Methods  for  Monotone  Nonlinear  Variational 
Inequalities,  Journal  of  Optimization  Theory  and  Applications,  Vol.  112,  pp.  111-128,  2002. 

13.  He,  B.S.,  Inexact  Implicit  Methods  for  Monotone  General  Variational  Inequalities,  Mathematical  Programming, 
Vol.  86,  pp.  199-217,  1999. 

14.  He,  B.S.,  A Class  of  Projection  and  Contraction  Methods  for  Monotone  Variational  Inequalities,  Applied 
Mathematics  and  Optimization,  Vol.  35,  pp.  69-76,  1997. 

15.  He,  B.S.,  Some  Predictor-Corrector  Methods  for  Monotone  Variational  Inequalities,  Technical  Report  95-65, 
Faculty  of  Technical  Mathematics  and  Informatics,  TU  Delft,  Delft,  Netherlands,  1995. 

16.  Hopfield,  J.J.,  Neural  Networks  and  Physical  Systems  with  Emergent  Collective  Computational  Ability, 
Proceedings  of  the  National  Academy  of  Sciences  of  the  USA,  Vol.  79,  pp.  2554-2558,  1982. 

17.  Hopfield,  J.J.,  Neurons  with  Graded  Response  Have  Collective  Computational  Properties  Like  Those  of  Two- 
State  Neurons,  Proceedings  of  the  National  Academy  of  Sciences  of  the  USA,  Vol.  81,  pp.  3088-3092,  1984. 

18.  [18]  Hopfield,  J.  J.,  and  Tank,  D.W.,  Neural  Computation  of  Decisions  in  Optimization  Problems,  Biological 
Cybernetics,  Vol.  52,  pp.  141-152,  1985. 

19.  Liao,  L.  Z.,  Qi,  H.  D.,  and  Qi,  L.Q.,  Neurodynamical  Optimization,  Journal  of  Global  Optimization  ( to  appear). 

20.  Leung,  Y.,  Chen,  K.Z.,  Jiao,  Y.C.,  Gao,  X.  B.,  and  Leung, K.S.,  A New  Gradient-Based  Neural  Network  for 
Solving  Linear  and  Quadratic  Programming  Problems,  IEEE  Transactions  on  Neural  Networks,  Vol.  12,  pp. 
1074-1083,2001. 

21.  Luenberger,  D.G.,  Linear  and  Nonlinear  Programming,  2nd  Edition,  Addi-son-Wesley  Publishing  Company, 
Reading,  Massachusetts,  1984. 

22.  Rockafellar,  R.T.,  Convex  Analysis,  Princeton  University  Press,  Princeton,  New  Jersey,  1970. 

23.  Rodr'iguez  V'azquez,  A.,  DomTnguez  Castro,  R.,  Rueda,  A.,  Huertas,  J.  L.,  and  S 'anchcz  Sinencio,  E., 
Nonlinear  Switch-Capacitor  Neural  Networks  for  Optimization  Problems,  IEEE  Transactions  on  Circuits  and 
Systems,  Vol.  37,  pp.  384-398,  1990. 

24.  Slotine,  J.  J.  E.,  and  Li,  W.,  Applied  Nonlinear  Control,  Prentice  Hall,  Englewood  Cliffs,  New  Jersey,  1991. 

25.  Xia,  Y.,  A New  Neural  Network  for  Solving  Linear  Programming  Problems  and  Its  Applications,  IEEE 
Transactions  on  Neural  Networks,  Vol.  7,  pp.  525-529,  1996. 


Impact  Factor  (JCC):  4.6257 


NAAS  Rating:  3.80 


